Neutron matter from chiral effective field theory interactions 
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The neutron-matter equation of state constrains the properties of many physical systems over a 
wide density range and can be studied systematically using chiral effective field theory (EFT). In 
chiral EFT, all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading 
order (N 3 LO). We present details and additional results of the first complete N 3 LO calculation of 
the neutron-matter energy, which includes the subleading three-nucleon as well as the leading four- 
nucleon forces, and provides theoretical uncertainties. In addition, we discuss the impact of our 
results for astrophysics: for the supernova equation of state, the symmetry energy and its density 
derivative, and for the structure of neutron stars. Finally, we give a first estimate for the size of the 
N 3 LO many-body contributions to the energy of symmetric nuclear matter, which shows that their 
inclusion will be important in nuclear structure calculations. 
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I. INTRODUCTION 

Chiral effective field theory (EFT) provides a system- 
atic expansion for nuclear forces including theoretical un- 
certainties [1] , where the development and applications of 
three-nucleon (3N) forces are a frontier [2J. In this con- 
text, neutron matter constitutes a unique laboratory for 
chiral EFT, because all many-body forces are predicted 
to N 3 LO [3]. This offers the possibility to provide reliable 
constraints based on chiral EFT interactions for neutron- 
rich matter in astrophysics, for the equation of state, the 
symmetry energy and its density dependence, and for the 
structure of neutron stars |U|S], but also allows to test the 
chiral EFT power counting and the hierarchy of many- 
body forces over a wide density range. In addition, the 
prediction of many-body forces makes neutron-rich nu- 
clei very exciting to test chiral EFT interactions against 
experiments at rare isotope beam facilities (6H2T]. 

Neutron matter has been studied in chiral EFT us- 
ing lattice simulations |22| and based on in-medium chi- 
ral perturbation theory |231 121] . In addition, neutron 
matter has been calculated using renormalization-group- 
evolved chiral EFT interactions [4 , where the renormal- 
ization group (RG) evolution improves the convergence 
of the many-body expansion around the Hartree-Fock en- 
ergy [25, 26J; and in a chiral Fermi liquid approach |27| . 
These studies demonstrated that 3N forces are significant 
at nuclear densities and that the dominant uncertainty is 
due to the truncation of 3N forces at the next-to-next-to- 
leading-order (N 2 LO) level [I]. Moreover, first Quantum 
Monte Carlo calculations with chiral EFT interactions 
are providing nonperturbative benchmarks for neutron 
matter at nuclear densities 128 . 
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Motivated by these studies and by the derivation of 
the parameter-free N 3 LO 3N and four-nucleon (4N) in- 
teractions [2TJ1 - |3"5] . we recently presented the first calcu- 
lation of the neutron-matter energy that includes all two- 
nucleon (NN), 3N and 4N forces consistently to N 3 LO 0. 
In this paper, we discuss details of our complete N 3 LO 
calculation and present additional results as well as appli- 
cations to astrophysics, for the equation of state and for 
the mass-radius relation of neutron stars. In addition, we 
give a first estimate for the size of the N 3 LO many-body 
contributions to the energy of symmetric nuclear matter 
in the Hartree-Fock approximation. This only presents a 
first step towards a complete calculation of nuclear mat- 
ter, where contributions from many-body forces beyond 
Hartree Fock are considerably more important than for 
neutron matter |34| . Our first results show that the in- 
clusion of N 3 LO 3N forces will be important in nuclear 
structure calculations. 

This paper is organized as follows. In Section [TT| we 
discuss the chiral EFT interactions included in this work. 
Details of the many-body calculation and convergence are 
given in Section |III| Our results for neutron matter are 
presented in Section |IV| including a detailed discussion 
of the uncertainties. In Section [V] we apply our results 
to the equation of state, in particular to the symmetry 
energy and its density dependence, and discuss the re- 
sulting constraints for the structure of neutron stars. We 
show first results for the N 3 LO 3N and 4N contributions 
in symmetric nuclear matter at the Hartree-Fock level in 



Section VI Finally, we summarize and give an outlook. 
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II. CHIRAL EFT INTERACTIONS 
A. N 2 LO and N 3 LO NN forces 

The largest interaction contributions to the neutron- 
matter energy arise from NN forces. For our past ap- 
plications of chiral EFT interactions to nucleonic mat- 



TABLE I. Spin-independent and spin-dependent two-body 
contact couplings Cs and Ct, respectively, for the N 3 LO NN 
potentials of Refs. [36 38J. 



TABLE II. Values of couplings ci, c 3 for the different NN 
potentials, as well as from Krebs, Gasparyan, and Epelbaum 
(KGE, Ref. [43]), and the range adopted in this work. 



NN potential 


C s [fm 2 ] 
-4.19 


C T [fm 2 ] 
-0.45 




ci [GcV- 1 ] 


c 3 [GeV" 1 ] 


EGM 450/500 MeV [36] 


N 2 LO/N 3 LO EGM NN [35] [36] 


-0.81 


-3.40 


EGM 450/700 MeV [3§] 


-4.71 


-0.24 


N 3 LO EMNN [33I3H] 


-0.81 


-3.20 


EM 500 MeV [33131] 


-3.90 


0.22 


N 2 LO KGE g3] 


-(0.26-0.58) 


-(2.80-3.14) 


EGM 550/600 MeV [35] 


-1.24 


0.36 


'N 2 LO' KGE (recom.) g3] 


-(0.37-0.73) 


-(2.71-3.38) 


EGM 600/600 MeV [36] 


3.45 


2.07 


N 3 LO KGE [43] 


-(0.75-1.13) 


-(4.77-5.51) 


EGM 600/700 MeV [36] 


1.31 


1.00 


N 2 LO this work 


-(0.37-0.81) 


-(2.71-3.40) 


EM 600 MeV [33 [35] 


-3.88 


0.28 


N 3 LO this work 


-(0.75-1.13) 


-(4.77-5.51) 



ter [H [31], the RG evolution has been used to evolve 
NN potentials to low-momentum interactions to improve 
the many-body convergence |25l 125] . In this work, we 
present calculations based directly on chiral EFT inter- 
actions without RG evolution and study the perturbative 
convergence following Ref. [3J. 

We investigate all existing NN potentials at N 2 LO and 
at N 3 LO of Epelbaum, Glockle, and Meifiner (EGM) [35J 
[35] with cutoffs A/A = 450/500, 450/700, 550/600, 
600/600 and 600/700 MeV, where A and A denote the 
cutoff in the Lippmann-Schwinger equation and in the 
two-pion-exchange spectral-function regularization, re- 
spectively; as well as the available N 3 LO NN potentials of 
Entem and Machleidt (EM) [371 EH] with cutoffs A = 500 
and 600 MeV. The EM 500 MeV potential is most com- 
monly used in nuclear structure calculations, while the 
EGM potentials have only been studied in some many- 
body calculations [33], although they allow to explore a 
wider cutoff range. 

The N 3 LO 3N and 4N forces involve the momentum- 
independent NN contact interactions Cs + Ct<Ti • 02 • bi 
particular, they mainly depend on Ct- The Cs, Ct val- 
ues of the different N 3 LO NN potentials are listed in Ta- 
ble II] for the neutron-proton case (the charge dependence 
contributes to higher-order charge-dependent 3N forces) . 
For a perturbative calculation, we require Wigner sym- 
metry (Ct = 0) to be fulfilled approximately at the in- 
teraction level. This is not the case for the EGM poten- 
tials with cutoffs 600/600 and 600/700 MeV, which have 
large spin-dependent couplings Ct ~ Cs (and even a re- 
pulsive spin-independent Cs), and would lead to large 
C T -dependent 3N forces at N 3 LO. 



B. N 2 LO 3N forces 

Three-nucleon forces enter at N 2 LO in the chiral EFT 
expansion without explicit Deltas [351 14"0"] . Due to the 
Pauli principle and the coupling of pions to spin, only the 
c\ and C3 parts of the long-range two-pion-exchange 3N 
interactions contribute at N 2 LO [4] (see also Ref. |41|). 
The same c, couplings also enter NN interactions at 
N~LO and have been determined from pion-nucleon or 



NN scattering. The c\, C3 values used in chiral NN po- 
tentials are given in Table [LT] Note however that the range 
adopted in the NN potentials of Table [IT] does not reflect 
the allowed range for the ci couplings, which are not sat- 
isfactorily constrained at present, e.g., with a range of 
c 3 = —(3.2-5.9) GeV -1 from different theoretical analy- 
ses (see Table I in Ref. [2]). 

We see from Table |LT| that, while the c\ value is of nat- 
ural size, the C3 value is large. This is due to the single- A 
excitation, which enhances C3 ~ 1/(itia — fn) by the A- 
nucleon mass difference to a large value (c\ — for a 
single- A excitation). In chiral EFT with explicit Deltas, 
the single-A contribution would in fact be included at 
one order lower; at next-to-leading order (NLO) in this 
case. The large C3 value has two effects. First, it leads 
to a slower convergence at the order when the c% con- 
tributions enter. This corresponds to topologies where 
A excitations are important. This can already be seen 
in the convergence pattern with NN interactions, where 
the leading two-pion-exchange NN interaction at NLO 
receives large contributions due to the large C, that en- 
ter the subleading two-pion-exchange NN interaction at 
N 2 LO [TJS2]. Therefore, for 3N and 4N forces important 
contributions to the N 3 LO interactions studied here can 
be expected in topologies where the c 2 ; couplings enter 
at N 4 LO |43l H4"] . This convergence pattern can be im- 
proved by including the Delta explicitly in chiral EFT. 
Second, the large C3 coupling in the N 2 LO 3N interaction 
also worsens the perturbative convergence of the many- 
body expansion around the Hartree-Fock energy. This is 
most important for the large C3 values considered in the 
N 3 LO calculation of this work (see Table [nj) . 

In addition, we list in Table [TTJ the Ci values extracted 
from a high-order analysis up to N 4 LO of Krebs, Gas- 
paryan, and Epelbaum (KGE, Ref. [13]). The KGE 
ranges at N 2 LO and N 3 LO are given in Table [jlj in ad- 
dition to values recommended to be used in an N 2 LO 
calculation that are tuned to capture the higher-order 
result. In this work, we take the KGE recommended 
Ci range for the N 2 LO calculation, minimally enlarged 
to include the Cj values of the NN potentials, and the 
KGE N 3 LO Ci range for our complete N 3 LO calculation. 
Note the large C3 value for the latter, which is still in the 
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FIG. 1. Diagrams contributing to the Hartree-Fock energy. 
These include the kinetic energy _Ekin and the first-order NN, 
3N, and 4N interaction energies E 
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range of Table I in Ref. |2|. We thus explore Cj values 
in the many-body interactions without varying the a in 
the NN potential. This is because changing the Cj in the 
NN potential would also require an adjustment of other 
couplings in the fit to NN data. We expect that some 
of the changes can be absorbed by the N 3 LO NN con- 
tact interactions, but it is very important to develop new 
N 3 LO NN potentials that can explore this sensitivity. 



C. N 3 LO 3N and 4N forces 

The many-body forces at N 3 LO are predicted by cou- 
plings in previous orders of the chiral EFT expansion. 
Hence, there are no new parameters for N 3 LO 3N and 
4N interactions [TJ. The subleading N 3 LO 3N forces have 
been derived recently |29ff3"T] . They can be grouped into 
five topologies, where the latter two depend on the NN 
contact couplings Ct and Cs (see the Appendix): 



T/ N J LO _ T/27T 

V 3N — V 



V 1 



+ yring + yl 



+ V 1/m . (1) 



Here, V 2 * ', V 2 ^' 1 ^ , and V nng denote the long-range 
two-pion-exchange, the two-pion-one-pion-exchange, and 
the pion-ring 3N interactions, respectively [30J. The 
terms y 27r - cont and V x l m are the short-range two-pion- 
exchange-contact 3N interaction and 3N relativistic cor- 
rections, respectively [21] • The latter are small [5] and 
depend also on the constants /3§ and J3g , which need to be 
chosen consistently with the unitary transformation used 
for the NN potentials [31] • In addition, there could be 
short-range one-pion-exchange-contact 3N interactions, 
but they have been shown to vanish at N 3 LO |31| . 

According to the chiral power counting, 4N forces en- 
ter at N 3 LO. They have been derived in Refs. [3J] [33J 
and depend also on the contact coupling Ct, but in neu- 
tron matter the CV-dependent parts do not contribute. 
There are seven 4N topologies that lead to non- vanishing 
contributions. In neutron matter only two three-pion- 
cxchangc diagrams (in Ref. [32 named V a and V e ) and 
the pion-pion-interaction diagram (V*) contribute |3J. 



III. MANY-BODY DETAILS 
A. Hartree Fock 

We calculate the energy per particle at the Hartree- 
Fock level and include contributions beyond Hartree Fock 



using many-body perturbation theory |3] [23] I34j . The 

Hamiltonian is given by H = T + Vnn + V3N + VW, 
where T is the kinetic energy and Vnn > VJjn > an d I4n de- 
note the NN, 3N, and 4N interactions, respectively. The 
Hartree-Fock contributions are shown diagrammatically 
in Fig. [I] At this level, the contribution of the A-nucleon 
interaction to the energy per particle is given by 
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with density n and short-hand notation i = kjO^. Here, 
Aa denotes the A-body antisymmetrizer and n^ — 
0(kp — ki) the Fermi-Dirac distribution at zero tem- 
perature. For the many-body forces, we use a Jacobi- 
momenta regulator. In terms of k^, this is given by 



f R _ e -p?+...+fci-k 1 -k 2 -...-k A _ 1 -k A )/(AA 2 )]"'= 



(3) 



where we take n oxp = 4 and consider 3N/4N cutoffs A = 
2 — 2.5 fm -1 . For the evaluation of 3N/4N forces, we 
use for the nucleon and pion mass, m — 938.92 MeV and 
7^ = 138.04 MeV, for the axial coupling g A — 1.29, and 
for the pion decay constant /„. = 92.4 MeV |30H33l 140] . 

As an example, we present details of the derivation 
of the Hartree-Fock energy from the N 3 LO two-pion- 
exchange 3N interactions. Their contributions can be 
grouped into two parts: one that shifts the c, couplings 
of the N 2 LO 3N forces and a part 



x/(4) _ 9a V^ 

v 2tt — 
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2567r/6^ fc (qf + m2)(q2 + m 2) 
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m v {ml + 3q 2 + 3q 2 + 4q, ; • q,-) 


+ (2ml + q?+q]+2q i -q J ) 


x (3ml + 3q 2 + 3q 2 + 4q, • C^)A(q k ) 





For the isospin part we have used that for neutrons 



(nnn\ Ti ■ Tj \nnn) = 1 



(nnn\ Ti ■ Tj x Tu \nnn) = , 



(4) 



(5) 
(6) 



and introduced the function F% w (qj,qi)> which absorbs 
all parts of the interaction except for the spin depen- 
dences. Furthermore, q^ = k^ — k^ and for F%J we use 
Qi + Qa = — *42 due to momentum conservation. Since 
the particles i, j, k are all neutrons and we sum over all 
possible spin states, the six different terms in the sum 
lead to identical contributions and we can write 



V, 



(4) 



2;r 



6(<Ti -qi)(o-3 •q 3 )F 2 ^ ) (qi,q 3 ) 



(7) 



For the spin trace Tr CT (123\A 3 V^ |123) we use that 
Pauli matrices are traceless and the relation uf(J b = 
S +ie a af. Thus, only the parts of the antisymmetrizer 
that contain the same-particle Pauli matrices as the po- 
tential need to be considered. In this case, the terms must 
contain or\ and <x 3 but not cr 2 . The antisymmetrizer is 
given by 



-4:3 



1 — p 12 — P 13 — p 23 -|- P 12 P 23 + P\ 3 P: 



23 : 



(8) 



with Pij = Ph — *^'' <T 3 ; where Pk exchanges the mo- 
menta of particles i and j. The last two terms can be 
written as 



P12P23 = 7^4(1 + tri ■ CT 2 + CT 2 ■ <T 3 



+ ct\ ■ cr 3 + icr 1 ■ cr 3 x <x 2 ) , 

^13^23 = ^iV^C 1 + (Tl ■ (T 2 + (T 2 ■ (T 3 

+ (Ti ■ ct 3 + icr 1 ■ cr 2 X cr 3 ) . (9) 
Thus, the only relevant terms of the antisymmetrizer are 



pk 

■M3 



pk pk 
^12^23 



pk pk 
■ r 13- r 23 



CT 1 ■ U 3 



(10) 



Multiplying this spin part with the potential leads to 



Tr„ 



ti • * 3 V^ 



= Tr 
= Tr„ 



6(5 ab + ie abd af)(S ac + ie ace a e 3 ) 



g^F^q^q,) 



(11) 



All terms containing Pauli matrices vanish when taking 
the trace, so that 



(4) _ 



(4) 



Tr CT <ri><r 3 V£> = 8 • 6q x • q 3 F™ ( qi) q 3 ) . (12) 



Thus, we obtain 



pk pk pk pk pk 

T, ,..4 3 vi 4) = 8 • 6 ( -tf + ^3 + ^3 



(4), 



x qi -qaF^'iqx,^) 



(13) 



Putting everything together yields for the spin-summed 
antisymmetrized matrix element 



(Vl 4) ) = I Tr ff (123| A 3 V 2 ^ |123) = 8 (123| (-& + **£& 



r>k r>k 
^13^23 



qi -q3F 2 (4) (q 1)q 3)|123) 



= -4 k 31 • kia F^ (kai , k 13 ) + 2 k 21 • k 13 F 2 (4) (k 21 , k 13 ) + 2 k 31 • k 23 F%? (k 31 , k 23 ) , 



P(4) 



= 4 ^ 3 F 2 ^(-k 13 ,k 13 )-k 12 -k 13 F 2 ^(-k 12 ,k 13 ) 



(4), 



(14) 



where k^ = k^ — kj, and we have relabeled the momen- 
tum indices in the last step, because the momentum in- 
tegrals are equal for the three neutrons and the regulator 
is symmetric under exchange of the momenta. 

Analogously, we obtain the expressions for the other 
N 3 LO 3N- and 4N-interaction matrix elements at the 
Hartree-Fock level. They are given in Appendix [X] The 
analytic derivations have been checked independently 
and using an automated Mathematica routine for the 
spin traces. 



B. Beyond Hartree Fock 

For nucleonic matter based on chiral EFT interactions, 
contributions beyond the Hartree-Fock level are impor- 
tant [H |23J Eii] . The dominant contribution to the energy 
is due to NN-NN correlations (E\ ). In addition, there 



(2) 



are NN-3N correlations (Ef> and E\>), 3N-3N correla 



r&h 



tions [E A 



(2) 



and E\ '), where the E, 



(2) 



follow the notation 



of Fig. |2j as well as NN-4N, 3N-4N and 4N-4N correla- 



tions. Based on the results of Refs. pUH], we expect the 

(2) 
residual 3N-3N contribution E§ and all contributions 

including 4N interactions to be small. 

The second-order contribution to the energy due to 
NN interactions and including 3N interactions as density- 
dependent two-body interactions is given by 



xx^ifnE/ 



d 3 k t 

(2tt) e 



:i2|K ( s 2) |34) 
-™k 4 ) 



1 a 

^ki"k 2 (l -n k3 )(l 

X 

£ki + £\a 2 ~ £ k 3 — £k 4 

x (2 7 r) 3 ,5(k 1 + k 2 - k 3 - k 4 ) . 



(15) 



.(2) 



where V as = (1 — -Pi2)Vnn + V3N * s * ne anti-symmetrized 
two-body interaction, which includes NN interactions 
and density-dependent two-body interactions from N 2 LO 
3N forces [¥. The latter are obtained by summing the 




FIG. 2. Second-order contributions to the energy due to NN- 



,-(2) 



n(2) 



NN correlations, E[ ' , NN-3N and 3N-3N correlations E, /3 

and E 4 , where the 3N forces enter as density-dependent two- 
body interactions, as well as the residual 3N-3N contribution 
Ei. not considered here. 



third particle over the occupied states in the Fermi sea 



V 



3N 
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C3 



d3k 3 „ a t/N^LO 

j2nf n ^ A3V ™ 



(16) 



At third order, we include particle-particle diagrams as 
in Ref . [34 . Their size provides a test of the convergence 
of the many-body calculation. We divide the third-order 

(3) 

particle-particle contributions into classes E\ , which are 

based on the E\ of Fig. ^ by adding one additional 
ladder and vertex with anti-symmetrized effective two- 



body interactions Vas ; = (1 — Pi2)Vnn 

(2) 

different diagrams E\ ' . 



^3N to the 



C. Convergence 

To study the perturbative convergence of the differ- 
ent NN potentials, we calculate the Hartree-Fock as well 
as second- and third-order energies with both free and 
Hartree-Fock single-particle energies. First, we consider 
NN interactions only and then study the changes when 
including also N 2 LO 3N forces. The results are shown 
in Figs. [3] and [4] respectively. The bands at each order 
range from using a free to a Hartree-Fock single-particle 
spectrum. In addition, we give in Table |III| the max- 
imal difference between the Hartree-Fock-spectrum re- 
sults at second order and those with a free or Hartee- 
Fock spectrum at third order for nuclear saturation den- 
sity no = 0.17 fm _ (corresponding to a Fermi momen- 
tum kp = 1.7 fm~ ). We take this energy difference as a 
measure of convergence for the potentials, as it includes 
both the uncertainty due to different single-particle ener- 
gies as well as the uncertainty in the convergence of the 
many-body calculation. 

At the NN level in Fig. |J the N 3 LO EGM potentials 
with cutoffs 450/500, 450/700, 550/600 and 600/600 MeV 
and the N 3 LO EM 500 MeV potential exhibit only small 
energy changes from second to third order. The larger- 
cutoff potentials (N 3 LO EGM 600/700 MeV and N 3 LO 
EM 600 MeV), however, show large changes from second 



TABLE III. Maximal energy difference between the second- 
and third-order contributions using a Hartree-Fock spectrum 
for the second-order and a free or Hartree-Fock spectrum 
for the third-order calculation at saturation density. Results 
are given for the different N 3 LO NN potentials at the NN- 
only level and including the leading N 2 LO 3N forces with 
A = 2.0 fin -1 , ci = -0.75 GeV" 1 and c 3 = -4.77 GeV -1 . 
The first three potentials exhibit a good convergence pattern 
with both NN-only and including N 2 LO 3N forces, and are 
therefore included in our complete N 3 LO calculation. Note 
that the N 3 LO EGM 600/600 and 600/700 potentials will not 
be considered in our complete N 3 LO calculation, because they 
have large Ct couplings (see the discussion of Table III. 



N 3 LO NN potential 


\AE {2/3) 

1 "-^NN-only 1 


| A£ (2/3) | 


EGM 450/500 MeV 


0.8 MeV 


0.6 MeV 


EGM 450/700 MeV 


0.4 MeV 


0.4 MeV 


EM 500 MeV 


1.1 MeV 


1.7 MeV 


EGM 550/600 MeV 


1.0 MeV 


3.1 MeV 


EGM 600/600 MeV 


0.2 MeV 


1.5 MeV 


EGM 600/700 MeV 


11.4 MeV 


16.1 MeV 


EM 600 MeV 


7.7 MeV 


9.1 MeV 



to third order, as well as a large band for the range of 
single-particle energies (especially for the EM 600 MeV 
potential). This demonstrates that these potentials are 
nonperturbative, see also Table [ITT] 

The convergence pattern is similar when the lead- 
ing N 2 LO 3N forces are included. We show the re- 
sults at this N 3 LO NN and N 2 LO 3N level in Fig. || 
for a 3N cutoff A = 2.0 fin - and a particular choice of 
ci = -0.75 GeV -1 and c 3 = -4.77GeV~\ although the 
general picture is unchanged for other coupling values. 

We find almost no change in the convergence pattern 
of the N 3 LO EGM 450/500 and 450/700 MeV poten- 
tials; see Table |III| This indicates that these potentials 
are perturbative for neutron matter. For the N 3 LO EGM 
450/500 MeV potential, this is expected already from the 
small Weinberg eigenvalues in Ref. J2H|, which are a nec- 
essary condition for the perturbative convergence. The 
perturbative convergence is a result of effective-range ef- 
fects [25], which weaken NN interactions at higher mo- 
menta, combined with weaker tensor forces among neu- 
trons, and with limited phase space at finite density due 
to Pauli blocking [25]. For the EM 500 MeV potential 
the inclusion of the N 2 LO 3N forces decreases the uncer- 
tainty estimate from the different single-particle energies, 
but increases the difference between second and third or- 
der. This can be seen comparing Figs. [3] and [4] and is 
reflected in the uncertainty estimate given in Table |III| 
Since this potential is most commonly used in nuclear 
structure calculations, we have decided to keep it in our 
complete N 3 LO calculation, in addition to the lower cut- 
off N 3 LO EGM 450/500 and 450/700 MeV potentials. 

The N 3 LO EGM 550/600 MeV potential is not used 
in the following calculations because its uncertainty esti- 



mate (see Table III I increases by a factor of three when 
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FIG. 3. (Color online) Energy per particle as a function of density for the different N 3 LO NN potentials of Refs. [35H38J . 
The dashed lines are Hartree-Fock results only. The filled and shaded bands are second- and third-order energies, respectively, 
where at each order the band ranges from using a free to a Hartree-Fock spectrum. 
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FIG. 4. (Color online) Energy per particle as a function of density for the different N 3 LO NN potentials of Refs. [35-35] and 
including the leading N 2 LO 3N forces. The dashed lines are Hartree-Fock results only. The filled and shaded bands are second- 
and third-order energies, respectively, where at each order the band ranges from using a free to a Hartree-Fock spectrum. All 
calculations are performed for a 3N cutoff A = 2.0 fm -1 and low-energy couplings c\ = 0.75 GeV -1 and cs = 4.77 GeV -1 . 



the N 2 LO 3N forces are included. This leads to a worse 
convergence pattern compared to the low-cutoff EGM po- 
tentials. For the N 3 LO EGM 600/700 MeV and EM 600 
MeV potentials we find the situation unchanged when in- 
cluding 3N forces, and thus, do not use these potentials 
for the following calculations. Even though the N 3 LO 
EGM 600/600 MeV potential exhibits a good convergence 
pattern, we will not use this interaction because it breaks 
Wigner symmetry at the interaction level (see the discus- 
sion of Table II]). Finally, we note that our findings for the 
EM 500 MeV potential are consistent with Ref. [JS] (see 
Fig. 6 therein), where the authors studied this potential 
at third order employing a Hartree-Fock spectrum. 



IV. RESULTS AND DISCUSSION 

Next, we present results using the EGM potentials with 
cutoffs 450/500 and 450/700 MeV and the EM 500 MeV 
potential. We discuss the individual contributions first 
and then show the complete N 3 LO results. 



A. N J LO NN and N 2 LO 3N forces 

The N 3 LO NN and N 2 LO 3N forces have been evalu- 
ated at the Hartree-Fock level and including second- and 
third-order contributions. Beyond Hartree Fock, N 2 LO 
3N forces are taken into account as density-dependent 
two-body interactions |3]. The kinetic energy, Hartree- 
Fock and individual higher-order interaction contribu- 
tions for the N 3 LO NN and N 2 LO 3N parts are given 
in Table [TV| for different values of A and the a couplings. 
Vanishing a in the 3N forces correspond to NN forces 
only. Table IV shows that the dominant higher-order 



contributions are due to the second-order NN-NN part 
E[ 2) . The second-order NN-3N parts E { 2 2) + E {2) are of 
the order of 1 MeV and only larger for the large 3N cutoff. 
All higher-order contributions with 3N forces are system- 
atically smaller. We emphasize that the N 2 LO 3N contri- 
butions beyond Hartree Fock are larger than in Ref. [4J, 
and therefore also the many-body calculation converges 
slower, because N 2 LO 3N forces are stronger due to the 
large N 3 LO values of the Cj couplings. 

The NN-only energies per particle are 14.7, 12.1, and 
12.9 MeV at saturation density for the EGM 450/500, 
450/700 MeV, and EM 500 MeV N 3 LO potentials, re- 
spectively. Inclusion of 3N forces at N 2 LO adds another 
7± 1.5 MeV per particle at saturation density (using the 
larger N 3 LO c, values, see Table |ll| . 



B. N 3 LO 3N and 4N forces 

The N 3 LO many-body forces have been evaluated in 
the Hartree-Fock approximation. We have not calcu- 
lated higher-order contributions because of their involved 
structure. The Hartree-Fock approximation is expected 



to be reliable based on the findings of Ref. [1] . In addi- 
tion, higher-order contributions with N 3 LO many-body 
forces are not enhanced by large c, couplings, and the 
N 3 LO many-body forces are smaller than at N 2 LO, lead- 
ing to smaller higher-order corrections. 

We show the individual contributions of the 3N and 
4N forces in Fig. [5] The bands correspond to the cutoff 
variation A = 2 — 2.5 fm~ . In the shorter-range two- 
pion-exchange-contact and the relativistic corrections 3N 
forces, three different bands are shown. These correspond 
to the different NN contacts, Ct and Cs, determined 
consistently for the different N 3 LO EM/EGM potentials. 

The two-pion-exchange 3N forces at N 3 LO yield an 
energy per particle of —1.5 MeV at saturation density, 
which is ~ 1/3 of the 3N contributions at N 2 LO and 
sets the natural scale. The two-pion-one-pion-exchange 
and the pion-ring 3N forces lead to relatively large con- 
tributions of —3.5 MeV and +3.3 MeV per particle at 
no, respectively. The contributions of the two-pion- 
exchange-contact 3N forces range between —2.8 MeV and 
+1.3 MeV per particle at no, depending on the NN poten- 
tial. In the topologies with relatively large expectation 
values, the large Ci couplings will enter in many-body 
forces at N 4 LO [43J. This may reflect important A con- 
tributions shifted to N 4 LO, as discussed above. Finally, 
the relativistic corrections contribute —(0.1 — 0.3) MeV 
to the energy per particle at no and are small compared 
with the other topologies. 

As shown in Fig. K3] (second panel) the sum of the N 3 LO 
3N contributions yields an energy of —(3 — 5) MeV per 
particle at saturation density for the EGM potentials and 
a small contribution of —0.5 MeV for the EM potential. 
This shows that the N 3 LO 3N contribution can be signif- 
icant, compared to the N 2 LO 3N energy of 7 ± 1.5 MeV 
per particle (note that the first panel of Fig. [6] only gives 
this contribution at the Hartree-Fock level). The rela- 
tively large N 3 LO 3N contributions are compensated by 
the larger N 3 LO Ci values, entering the 3N force at N 2 LO. 
This can be seen in Fig. [6] where the total 3N contribu- 
tion at N 3 LO (third panel) is compared at the Hartree- 
Fock level to the 3N contribution at N 2 LO (fourth panel), 
which uses the Cj values recommended for an N 2 LO cal- 
culation (see Table [h| . For the EGM potentials the total 
3N contribution changes by less than 1 MeV going from 
N 2 LO to N 3 LO. Because the N 3 LO 3N contribution is 
small for the EM potential, this results in a difference of 
about 3 MeV when going from N 2 LO to N 3 LO for the 
EM case, due to the modified c, couplings at N 3 LO. 

Only three N 3 LO 4N topologies give nonvanishing 
contributions to neutron matter. We show their re- 
sults in Fig. [5j The two three-pion-exchange diagrams 
V a and V e are attractive with energies of —0.16 MeV 
and —0.25 MeV per particle at saturation density. The 
pion-pion-interaction 4N forces (V*) are repulsive with 
0.22 MeV per particle at no. The latter two diagrams al- 
most cancel each other, such that the total contribution 
of the leading 4N forces is about —0.18 MeV per particle 
at no- However, also for the 4N forces additional larger 



TABLE IV. Contributions from different N 3 LO NN potentials and the leading N 2 LO 3N 
per particle in MeV at nuclear saturation density. A Hartree-Fock spectrum for the single- 
3N force is for different A in fm -1 and for different ci, C3 in GeV -1 . 
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FIG. 5. (Color online) Energy per particle as a function of density for all individual N 3 LO 3N- and 4N-force contributions 
to neutron matter at the Hartree-Fock level. The bands are obtained by varying the 3N/4N cutoffs A = 2 — 2.5 fm" 1 . For 
the two-pion-exchange-contact and the relativistic-corrections 3N forces, the different bands correspond to the different NN 
contacts, Ct and Cs, determined consistently for the N 3 LO EM/EGM potentials. The inset diagram illustrates the 3N/4N 
force topology of the particular contributions. 
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FIG. 6. (Color online) Contributions from 3N forces in the Hartree-Fock approximation at N 2 LO plus N 3 LO (first three panels) 
in comparison with the 3N contribution in a N 2 LO calculation (fourth panel). The first panel shows the N 2 LO 3N contribution 
in the N 3 LO calculation, using the N 3 LO values of the a couplings, and the second panel gives the N 3 LO 3N contribution. 
The third panel shows the total 3N contribution at N 3 LO (the sum of the first two panels). This is compared in the fourth 
panel to the 3N contribution at N 2 LO, using the d values recommended for an N 2 LO calculation (see Table lllk. For the EGM 
potentials the total 3N contribution at N 3 LO differs by less than 1 MeV compared to the N 2 LO results. However, for the EM 
potential, the result changes by almost 3 MeV. All bands include the a range from Table |TT| and the 3N cutoff variation. 



contributions from A excitations may arise at N 4 LO [H] ■ 
At the Hartree-Fock level, the 3N/4N contributions 
change by less than 5% if the cutoff is taken to infinity 
(i.e., fji — 1). However, since we also include N 2 LO 3N 
forces beyond Hartree Fock, a consistent regulator is re- 
quired. Finally, we compare our 4N results with those of 
Refs. |44[I47|. which considered only the 4N interactions 
V e and V? and found their sum to be about —11 keV per 
particle at no. This is in agreement with our results, if 
we take f R = 1 as in Refs. |44ll47] . 



C. Complete calculation at N 3 LO 

The complete N 3 LO result for neutron matter is shown 
in Fig. [7] which includes all many-body interactions to 
N 3 LO |Sf . At saturation density, we obtain for the energy 
per particle 



— (n ) = 14.1- 21.0 MeV. 



(17) 



This range is based on different NN potentials, a varia- 
tion of the couplings c\ = —(0.75 — 1.13) GeV~ , C3 = 
-(4.77-5.51) GeV -1 , and on the 3N/4N-cutoff variation 
A = 2 — 2.5 fm - . In addition, the uncertainty in the 
many-body calculation is included, as discussed above. 

As shown in Fig. [71 our results are consistent with pre- 
vious calculations based on RG-evolved NN interactions 
at N 3 LO and 3N interactions at N 2 LO (4]. These calcu- 
lations adopted a conservative Cj range, but are based on 
the EM 500 MeV NN potential only, which results in a 
narrower band compared to the N 3 LO band. In Ref. [3], 
we compared our results to calculations based on lattice 
EFT 1 22 J and Quantum Monte Carlo at low densities [48] , 
as well as to variational methods [49 and Auxiliary Field 
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FIG. 7. (Color online) Neutron-matter energy per particle 
as a function of density including NN, 3N, and 4N forces to 
N 3 LO. The three overlapping bands are labeled by the differ- 
ent NN potentials and include uncertainty estimates due to 
the many-body calculation, the low-energy d constants, and 
by varying the 3N/4N cutoffs (see text for details). For com- 
parison, we show the results for the RG-evolved NN EM 500 
MeV potential including only N 2 LO 3N forces from Ref. |lj- 



Diffusion Monte Carlo [50] based on phenomenological 
NN and 3N potentials, and found that they are also con- 
sistent with the N 3 LO band. However, the latter calcu- 
lations do not provide theoretical uncertainties. 

In Fig. m we compare the convergence from N 2 LO to 
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FIG. 8. (Color online) Neutron-matter energy per particle as 
a function of density at N 2 LO (upper blue band that extends 
to the dashed line) and N 3 LO (lower red band). The bands 
are based on the EGM NN potentials and include uncertainty 
estimates as in Fig. T7\ 



N 3 LO in the same calculational setup. For this compar- 
ison, we only consider the EGM potentials with cutoffs 
450/500 and 450/700 MeV, since no EM N 2 LO poten- 
tial is available. This leads to an N 3 LO energy range of 
14.1 - 18.4 MeV per particle at n . For the N 2 LO band 
in Fig. [8] we have estimated the theoretical uncertainties 
in the same way, and find an energy of 15.5 — 21.4 MeV 
per particle at no- The two bands overlap but the range 
of the band is only reduced by a factor of 2/3, which is 
larger than the 1/3 expected from the EFT power count- 
ing. We attribute this to A effects (as discussed above). 
This can be improved by including the Delta in chiral 
EFT explicitly or by going to N 4 LO @5]. 

Finally, it is important to construct NN potentials at 
N 2 LO and N 3 LO covering the range of the Ci values. At 
N 3 LO, we expect that the differences in the Ci can be ab- 
sorbed partly by Q 4 contact interactions in the fits to NN 
scattering. In addition, the many-body-calculation un- 
certainties can be reduced further by including the N 3 LO 
many-body forces beyond the Hartree-Fock level. 



APPLICATIONS 



TABLE V. Ranges for the symmetry energy 
derivative L at nuclear saturation density. 


S v and its density 


range 


Symmetry energy 
Density derivative 


S v (no) 
L(n ) 




28.9 
43.0 


- 34.9 MeV 
-66.6 MeV 



follow Ref. [52_ and take an expression that includes ki- 
netic energy plus interaction energy that is quadratic in 
the neutron excess 1 — 2x, where x is the proton fraction, 



e(n, x) = T 



■ x* + (1 — x)s (2n)3 

— [(2a — 4a.£,)a;(l — x) + aj~\n 
+ [(27? - 4»7i)a;(l - x) + r] L ] n^ 



(18) 



where n = n/n a and T = (37r 2 no/2) 2 / 3 /(2m) = 
36.84 MeV is the Fermi energy of symmetric nuclear mat- 
ter at saturation density. The parameters a = 5.87 and 
rj = 3.81 are determined through fits to the empirical sat- 
uration point of nuclear matter, and a^ and t\l through 
fits to the neutron-matter results of Fi g. [7] (for details on 
this strategy, see Ref. [55]). Equation (18) provides very 
good fits to the N 3 LO energy band. 

We can then calculate the symmetry energy 



S v (n) 



1 d 2 e(n,x) 

8 



dx 2 



and its density derivative 



L(n) 



3 d 3 e(n, x) 
8 dndx 2 



-\,x=\ji 



n=l,x=l/2 



(19) 



(20) 



The L parameter basically determines the pressure of 



neutron matter. In addition, because the expression (18 1 



is fit to the empirical saturation point (with small uncer- 
tainties), the symmetry energy and its density derivative 
at n and their theoretical uncertainties are essentially 
determined by the neutron-matter results. 

The predicted ranges for S v and L at saturation density 
are given in Table W\ In Ref. (3] , we have shown that S v 
and L are also correlated and overlap with the results for 
RG-evolved NN interactions with N 2 LO 3N forces [3TJ 
IB"2"] . but due to the additional density dependences from 
N 3 LO many-body forces, this correlation is not as tight. 
The S v and L ranges are also in very good agreement 
with experimental constraints from nuclear masses |53| 
and from the dipole polarizability of 208 Pb [51] (see also 
Refs. [3]|5T]). 



A. Symmetry energy and its density derivative 

The symmetry energy S v and its density derivative L 
provide important input for astrophysics [51J. To calcu- 
late these, we need to extend the neutron-matter energy 
to asymmetric matter. For the energy per particle e, we 



B. Constraints for supernova equations of state 
and neutron stars 



The neutron-matter results also provide constraints for 
the nuclear equation of state. Here we focus on compar- 
isons to equations of state for core-collapse supernova 
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FIG. 9. (Color online) Comparison of the neutron-matter 
energy at N 3 LO of Fig. M (red band) with equations of state 
for core-collapse supernova simulations provided by Lattimer- 
Swesty (LS |55| with different incompressibilities 180, 220, 
and 375 MeV), G. Shen (FSU2.1, NL3 |BT]), Hempel (TM1, 
SFHo, SFHx [62j), and Typel (DD2 [56]). 



simulations. In Fig. |9| we compare the N 3 LO neutron- 
matter band (red band) to the Lattimer-Swesty (LS) 
equation of state |55| , which is most commonly used 
in simulations, and to different relativistic mean-field- 
theory equations of state based on the density functionals 
DD2 ISP, FSU2.1 [S7], NL3 |5S], SFHo, SFHx [55], and 
TM1 [60 . At low densities only the DD2, FSU2.1 and 
SFHx equations of state are consistent with the N 3 LO 
neutron-matter band. The other supernova equations 
of state underestimate the energy for densities below 
~ 0.5no and even at higher density in the LS cases. This 
density range covers the outer regions of the (proto-) neu- 
tron star, where also protons, nuclei, and electrons are 
relevant. Nevertheless, the deficiencies in the nuclear in- 
teractions of these equations of state will also affect the 
chemical potentials and the neutrino response. Around 
saturation density, the LS and SFHo equations of state 
become consistent with the N 3 LO band. We also find 
that the NL3 and TM1 equations of state have a too 
strong density dependence, which leads to unnaturally 
large S v and L values. In addition, Fig. [9] exhibits a 
strange density dependence of SFHx. 

Next, we use the N 3 LO neutron-matter results to pro- 
vide constraints for the structure of neutron stars. We 
follow Ref. [52J for incorportating beta equilibrium and 
for the extension to high densities using piecewise poly- 
tropes that are constrained by causality and by the re- 
quirement to support a 1.97±O.O4M neutron star |63| . 
the heaviest precisely measured neutron star to date. 



FIG. 10. (Color online) Constraints on the mass-radius di- 
agram of neutron stars based on our neutron-matter results 
at N 3 LO following Ref. |52j for the extension to neutron-star 
matter and to high densities (red band), in comparison to 
the constraints from calculations based on RG-evolved NN 
interactions (thick dashed blue lines) |52| . We also show the 
mass-radius relations obtained from the equations of state for 
core-collapse supernova simulations shown in Fig. M [551 1571 - 
|60]|6l[65]. The legend for the thin lines is as in Fig. [9] 



The resulting constraints on the neutron star mass-radius 
diagram are shown in Fig. [lO] by the red band. This 
band represents an envelope of a large number of individ- 
ual equations of state reflecting the uncertainties in the 
N 3 LO neutron-matter calculation and in the polytropic 
extensions to high densities [52] 
predicted radius range of Ref. [52] of 9.7 
a 1.4 M Q neutron star. The largest supported neutron 
star mass is found to be 3.1 M©, with a corresponding 
radius of about 14 km. We also find very good agree- 



Figure 10 confirms the 
13.9 km for 



ment with the mass-radius constraints from the neutron- 
matter calculations based on RG-evolved NN interactions 
with N 2 LO 3N forces [52], which are shown by the thick 
dashed blue lines in Fig. [10] 

In addition, we show in Fig. [10] the mass-radius rela- 
tions obtained from equations of state for core-collapse 
supernova simulations [55, 57 60] [64] [65]. The inconsis- 
tency in Fig. |9]of many of the equations of state with the 
N 3 LO neutron-matter band at low densities results in a 
large spread of very low mass/large radius neutron stars, 
where the red band is considerably narrower in Fig. [10] 
(note that the red band includes a standard crust equa- 
tion of state below 0.5 no [52]). For typical neutron stars, 
our calculations rule out the NL3 and TM1 equations of 
state, which produce too large radii. Finally, we em- 
phasize that these constraints are not only important for 
neutron star structure and for the supernova equation of 
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FIG. 11. (Color online) Energy per particle versus density for all individual N 3 LO 3N- and 4N-force contributions to symmetric 
nuclear matter at the Hartree-Fock level. The bands are obtained by varying the 3N/4N cutoff A = 2 — 2.5 fm -1 . For the 
two-pion-exchange-contact, the relativistic-corrections 3N forces, and the short-range 4N forces, the different bands correspond 
to the different NN contacts, Ct and Cs, determined consistently for the N 3 LO EM/EGM potentials. The inset diagram 
illustrates the 3N/4N force topology of the particular contribution. 



state, but also provide nuclear physics constraints for the 
gravitational wave signal in neutron star mergers [66II67J . 



VI. FIRST ESTIMATE FOR SYMMETRIC 
NUCLEAR MATTER 



We present first results for the N 3 LO many-body forces 
in symmetric nuclear matter in the Hartree-Fock approxi- 
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mation. However, we emphasize that these results should 
be considered as a preview and to show their importance, 
because it is crucial to include contributions beyond the 
Hartree-Fock level [34 . Such calculations can also be 
facilitated by a similarity RG evolution of NN and 3N 
forces |68[ 169] in order to improve the convergence of the 
many-body calculation. 

The energy per particle of symmetric matter is evalu- 
III A| summing also over both isospin 



Sect! 



ated a 

states [see Eq. (Bl )]. In Appendix |B| the expressions for 
the N 3 LO 3N- and 4N-interaction matrix elements are 
given in detail. Our results for the individual contribu- 
tions from N 3 LO many-body forces are shown in Fig. 11 



Compared to the neutron-matter results, the individual 
contributions are larger in magnitude in symmetric mat- 
ter, requiring calculations beyond the Hartree-Fock level. 
However, the bands from cutoff variation are narrower, 
because the Fermi momentum corresponding to satura- 
tion density is lower in symmetric matter. 

For the two-pion-exchange N 3 LO 3N forces the en- 
ergy is small with 0.24 MeV per particle at no, due to 
cancellations among the individual parts in symmetric 
matter. The other 3N topologies are large and attrac- 
tive: the two-pion-one-pion-exchange and the pion-ring 
3N interactions give energies of —6.5 MeV and —3.6 MeV 
per particle at no, respectively. The contribution of the 
two-pion-exchange-contact 3N interaction ranges from 
-7.0 MeV to +3.4 MeV, depending on the NN potential. 
As expected from our neutron-matter results, the large 
3N contributions in these topologies can be attributed to 
the physics from A excitations, which will lead to large 
Ci contributions at N 4 LO in these topologies (or at N 3 LO 
in A- full chiral EFT). As in neutron matter, the contri- 
butions from relativistic corrections 3N forces are small 
with -(0.24 - 0.39) MeV per particle at n . 

Since nuclear saturation is a result of cancellation ef- 
fects of large energy contributions [33], the increased 
strengths of the Cj couplings at N 3 LO compared to N 2 LO 
is expected to play an important role for predictions of 
symmetric matter. Furthermore, in contrast to neutron 
matter, we find that the total N 3 LO 3N contribution at 
the Hartree-Fock level depends more strongly on the NN 
potentials used: for the EM 500 MeV potential, we find 
—7 MeV per particle at no, whereas for the EGM poten- 
tials, we find —(15—17) MeV. To understand this better, 
improved NN potential fits (following Ref. [70]) and also 
for different c, couplings will be important. These N 3 LO 
energies should be compared with a total N 2 LO 3N en- 
ergy at the Hartree-Fock level of the order of 15 MeV per 
particle at n , using the large N 3 LO c* values (see Table |llj 
and accordingly chosen C4 = 3.34 — 3.71 GeV -1 [33]) and 
typical cd , ce values [2S] ■ All these findings point to that 
including N 3 LO 3N contributions beyond the Hartree- 
Fock level will be crucial. 



Figure 11 also shows our results for the individ- 
ual N 3 LO 4N-force contributions in symmetric matter. 
The long-range three-pion-exchange 4N interactions V a 
and V e are attractive with energies —0.32 MeV and 



—0.39 MeV per particle at no, respectively; while the V c 
interaction is repulsive with 0.21 MeV per particle at no- 
The pion-pion-interaction 4N force V* also gives a re- 
pulsive contribution of 0.33 MeV per particle at no- The 
shorter-range parts V k , V , and V n contain one or two 
spin-dependent NN contact interactions and depend on 
Ct- The two mid-range topologies involving only one 
NN contact (V k and V 1 ) almost cancel against each other 
(—0.11 to +0.05 MeV per particle at no, depending on the 
NN potential, and -0.05 to +0.10 MeV, respectively). 
The shortest-range topology with two NN contacts (V™) 
contributes even less (—0.06 to —0.01 MeV per particle 
at no). In total, the leading 4N forces give an attractive 
contribution of —(0.18 — 0.23) MeV per particle at no, 
with a strong density dependence ~ n 3 . 

As a check, we can compare our results to the studies of 
the V e and V f 4N forces of Refs. [33] SZ], which obtained 
a contribution to the energy per particle of —53 keV at no- 
This is in agreement with our result for the sum of these 
two topologies: —(56 + 2) keV, where the small difference 
is due to fn = 1 in Refs. |4"4l |4"7] . So far only the leading 
4N forces have been derived completely. Recently, Kaiser 
studied A contributions to 4N forces [33]) which enter at 
N 4 LO in A-less chiral EFT. Similar to the N 3 LO versus 
N 2 LO 3N forces, these contributions are enhanced by the 
large c, values, and Kaiser found for these partial N 4 LO 
4N contributions a larger energy of ~ 2 MeV per particle 
at saturation density. 

Finally, we compare our results for symmetric matter 
with first calculations of the 4N contributions to the 4 He 
ground-state energy. These were studied in Ref. [71] per- 
turbatively based on the same N 3 LO 4N forces. We agree 
with the sign of the 4N contributions for all topologies 
and obtain a similar total energy correction when taking 
a density ~ no/3. Also, the estimate of Ref. [72] for the 
V e 4N contribution to the 4 He ground-state energy gave 
— 18keV per particle, which is of the same order as our 
results at ~ no/3. 



VII. SUMMARY AND OUTLOOK 

We have presented details and additional results of the 
first complete N 3 LO calculation of the neutron-matter 
energy based on chiral EFT NN, 3N, and 4N interac- 
tions [3]- Our results for the energy per particle at sat- 
uration density give a range of 14.1 — 21 MeV, which in- 
cludes uncertainties from different NN potentials, from 
the C\ and C3 couplings in 3N forces (these dominate), 
from varying the cutoff in many-body forces, and from 
the uncertainties in the perturbative many-body expan- 
sion around Hartree Fock. For more systematic studies, 
it will be important to develop NN potentials that ex- 
plore the different a couplings. 

We have found large contributions to the energy from 
N 3 LO 3N forces in topologies where A excitations are im- 
portant. Therefore, an improved EFT convergence is ex- 
pected in chiral EFT with explicit A degrees of freedom. 
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In contrast, contributions from the leading 4N forces are 
found to be small (see also Refs. |44(l4"T]). We have pre- 
sented a first estimate for the N 3 LO many-body contri- 
butions to the energy of symmetric nuclear matter, where 
also large N 3 LO 3N forces and small leading 4N forces 
are found. Our results for symmetric matter show that 
the inclusion of N 3 LO 3N forces will be important in nu- 
clear structure calculations, and that it is crucial to go 
beyond the Hartree-Fock approximation. 

Recently, first Quantum Monte Carlo calculations with 
chiral EFT interactions are providing nonperturbative 
benchmarks for neutron matter and validate the pertur- 
bative expansion for chiral NN potentials with low cut- 
offs [28]. Extending these calculations to 3N forces and 
N 3 LO will be important. In addition, the many-body 
uncertainties can be reduced in the future by a similarity 
RG evolution of NN and 3N forces [Ml EH], which im- 
proves the many-body convergence and will also enable 
studies with the chiral NN interactions, which were found 
to be nonperturbative in the present calculations. 

In addition, we have discussed the impact of our re- 
sults for astrophysics: The predicted ranges for the sym- 
metry energy S v and its density derivative L are S v — 



28.9-34.9 MeV and L = 43.0-66.6 MeV, which are con- 
sistent with recent experimental constraints |51[ 153) 154] . 
Many of the equations of state for core-collapse supernova 
simulations were found to be inconsistent with the N 3 LO 
neutron-matter band. By extending our neutron-matter 
results to neutron-star matter and to high densities, we 
confirm the predicted radius range of 9.7 — 13.9 km for 
a 1.4 Mq neutron star |52) and find a maximal neutron 
star mass of 3.1 Mq. 
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Appendix A: N 3 LO neutron-matter matrix elements 



In this Appendix we present the 3N and 4N matrix elements defined as 



1 






cr lt ...,aA i-l=£...^i a 



l---A\A A J2 VA N (ii,...,i A )|l ■■■>!> , 



(Al) 



entering the neutron-matter Hartree-Fock calculation [see Eq. (pi)] of the N 3 LO many-body forces. 

We use the short-hand notation for the momentum transfer k^ = kj — kj, k(ij)(kl) = ky + kfc; and P^ = ' 2 J , 
and pion propagators K io = kfj + ml and K {ij){kl) = k 2 vj){kl) + m£. 



1. Two-pion-exchange 3N 



A 



W> = f [~2S Cl mi 



ki2 • k23 

K12K23 



k 2 

ft 12 

K 2 

-"-12 



<5cs 



(4) 



(ki2 • k 23 ) 2 kf 2 
K12K23 



k( 3 ^ i("ki3, k 13 ) - k 12 • k 13 F£ ^(-kia, k 18 ) , 



K 2 



with shifts in the low-energy couplings 6c% = — 0.13GcV and Sc 3 — 0.89 GeV 1 (see Ref. 



(A2) 
]) and the function 



^ii(qi,q 2 ) 



ig\ 



m n {ml + 3ql + Zql + 4qi • q 2 ) 



+ (2ml + Qi+ql + 2qi • q 2 )(3m2 + Z q 2 + 3g| + 4qj • q 2 )A{\q 1 + q 2 |) 



(A3) 



where A(q) = l/(2q) arcta,n[q/(2m 7V )] denotes the loop function 



2. Two-pion one-pion-exchange 3N 



15 



(v^ 1 ! 



Fi(h 



+F 4 (k 1 



(kis • k 13 ) 2 



K 13 



-"■13 -"23 



,( k 



12 ■ tVi.::r ,_. , , ^ k\ z 



K 



13 



+ F 5 {k 12 )-^ - F Q {k l2 ) 
A 13 



>-12 ' K-13 



K 



13 



f 3 (m|^-^(o)|^ 

A 13 -"23 



F 7 (fci 



^13 



A' 



13 



(A4) 



with structure functions Fi{q) to Ff(q) defined in Eqs. (2.17)— (2.20) of Ref. 



3. Pion-ring 3N 



(Ki lnSX 



'3N 



= 4 



-3 i?i(k 12 , 0) + 3 i?x(k 12 , k 23 ) - k\ 2 # 2 (k 12 , 0) + k\ 2 R 2 (k 12 , k 23 ) + k 12 • k 23 i? 3 (k 12 , k 23 ) 

k i2 • k 23 A 4 (ki 2 , k 23 ) + k\ 3 i? 5 (ki2, k 23 ) + 2A 6 (0, 0) - A 6 (k 12 , 0) - A 6 (0, k 12 ) - A 6 (-ki 2 , k 12 ) 
i? 6 (ki 2 , k 23 ) - k\ 2 i? 7 (-ki 2 , k i2 ) + k\ 2 R 7 (k 12 , k 23 ) + k\ 2 i? 8 (-ki 2 , k i2 ) + k i2 • k 23 A 8 (ki 2 , k 23 ) 
k\ 2 i? 9 (-k 12 , k 12 ) + k 12 ■ k 23 A 9 (k 12 , k 23 ) - 3i?i (-k 12 , k 12 ) + 3i? 10 (k 12 , k 23 ) + 25i(0, 0) 
Si(ki 2 , 0) - Si(0, k 12 ) - Si(-ki2,ki 2 ) + Si(ki2,k 23 ) - ^ 2 5 2 (-k 12 , k 12 ) + k 2 12 S 2 (k 12 , k 23 ) 
^^(-kia.k^) +k i2 • k 23 S' 3 (ki 2 ,k 23 ) + fci 2 S , 4 (-ki 2 ,ki 2 ) + k i2 • k 23 S , 4 (ki 2 ,k 23 ) 

^ 2 5 5 (-k 12 ,k 12 ) + fc 2 2 3 5 5 (k 12 ,k 23 ) - 35 6 (-ki2,k 12 ) + 3S 6 (ki2,k2 3 )l , (A5) 



where the structure functions Ri and Si are defined in Eqs. (A2) and (A7) of Ref. 



4. Two-pion-exchange contact 3N 



(% 



27r-cont\ 
3N / 



A. 



Cn 



A 






4to 2 



"-12 



~2{2ml + k 2 12 )A(k 12 ) 



^-(2ml + kl 2 )A(k 12 ) 



(A6) 



Relativistic-corrections 3N 



(V l/r ' 

\ K 3N 



k 2 2 Fl /m (k 12 ,k 1 2)+'ki2 ■ k 23 F} /m (k 12 ,k 23 ) - (ki 2 x k 23 ) 2 F 1 2 /m (-k 12 ,k 13 ,P 12 ,P 23 ) 
-k 2 12 F? /m (k l2 ,k 12 ) + k l2 ■k 2 3F 1 3 /m (ki 2 ,k 23 )- (ku x k 18 ) ■ (k xa x P 23 )^ 1 4 /m (k 12 ,k 13 ) 
- (ku x ki 3 ) • (ki 2 x P 13 )f 1 5 /m (k 12 ,ki 3 ) + fc? 2 Jf /m (ki 2 , k 23 ) - fc 2 2 ^ 7 /m (k 12 , -k 12 ) 
-k 2 12 F 7 1/m (k 12 ,k 23 ) - k 2 2 F? /m (k 12 ,P 12 ,P 23 ) - k 2 12 F? /m (k 12 ) + kl 2 FlJ m {k l2 ) + k\ 2 F\) m {k 12 ) 



(A7) 



with 
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pi , n]= A (l-2 / g 8 )(q 1 -q 2 ) 2 

1/mlqi ' q2j 16m/* (ql+mlYiql + ml) ' 

g\ (1 - 2/3 8 ) qi ■ q 4 + (1 + 2^ 8 ) qi ■ q 3 



^i 2 / m (qi,q2,q 3 ,q4) 



8m/4 ( 9 2 +m 2)2( g 2 +m 2) 



^i 3 /TO (qi,q2) = - 



gj (2fo - l)g 2 

16m/* (q 2 + m 2 )(<7 2 + m 2 ) 



<ii 



^i 6 /m (qi,q2) 



flA 



:C : 



(1 - 2^ 8 ) qi • q 2 



= -^i 4 / m (qi>q2)y = 

Cs 



- i? i 5 / m (qi.q2) 



g?(2 j 8 9 -l) 

2(2^ 9 + l) ' 



^i/ m (qi>q2)7^, 



4m/r 6 ' (^m 2 ) 2 '1/™^""^ 

p8 , x 9a n (1 ~ 2^ 8 )qi ■ q 3 + (1 + 2^ 8 ) qi ■ q 2 

^i/ m iqi>q2jq3j — — «^t 



m/2 



2^2 



A 



(qf + ml) 

Pl0 , x^S 



C, 



p 9 c^ - 9a r 2l3g ~ x - p 10 ivi^ 5 - p 11 (^ u ' 5 



(A8) 
(A9) 
(A10) 
(All) 
(A12) 
(A13) 



6. Three-pion-exchange and pion-interaction 4N 



<VTn> 






(kj x k 2 ) • k 34 + (k 3 x k 4 ) • k 12 ]" 



K ^Kf m23) K 24 K 12 K 2 4 K 34 K 12 K 2 3 K 14 



fc^ 4 (ki4 x k( 14 ) (23) ) 2 ki4 • k 24 (k( 14) ( 23 ) x k i4 ) • (k (14)(23 ) x k 24 ) 



K 2 K 2 
- n 'i4- fv (i 4 )( 23 ) 



K 14 K (U)(23) K2 -± 
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K 12 K 2 K 34 



K 12 K 2 K 14 



(V&) = 
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31 
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(A14) 



k 23 • k 24 k 23 • ki4 k i2 • k 24 
2 ki 3 • (ki 3 + k 24 ) + 2 k 34 • ki 3 + 2 k 34 • k 23 

-^13-^23-^24 ^14-^23-^34 -?^1 2 -^24-^34 
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K+2^ (1 2)(34))||| 2 



/ \ _ 9 A 

32/ 

-2(K U + A"(34)(21) + ^23 



31 
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14)(32 
12 J1 34 

k 42 • k 24 ki3 • k 34 
-fi'l 2 ^13-^24^34 



IK 



ki2 -k 34 ki 4 • k 



23 



13 J 



-^12-^14-^23-^34 



(A15) 
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Appendix B: N 3 LO symmetric nuclear-matter matrix elements 



We now turn to the 3N and 4N matrix elements defined as 



(^n) = J, Y, E £-MAa E V AN (h,...,i A )\l---A) 



Tl,...,TA V\, ■■-,<! A 



i 1 ^...^i / 



(Bl) 



entering the symmetric nuclear-matter Hartree-Fock calculation of the N 3 LO many-body forces. 



1. Two-pion-exchange 3N 
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(V, 



2tt\ 
3N/ 



P 

J IT 



Scirnl 



,2 c(4) 



ki2 ■ K.23 



ft 12 
L ^12^23 ' ~K\ 2 



Sc 3 
P 

J TV 

r(4) 



(kia • k 23 ) 2 
K\ 2 K 23 



-"-12 



<5c 4 (ki2 x k 2 



2k 2 3 F 2 ^ 1 (-k 13 ,k 13 ) - k 12 • ki3F£ xt-kia.ku,) - (k 12 x k 13 ) 2 F 2 ^ 2 (-k 12 ,k 13 



/ 2 #12^2 

^2 F (4) 



with shifts in the low-energy couplings dci, Sc 3 = — <5c 4 , the function F 2 J X is as given in Appendix |A| and 



F^(qi,q2) = - 



9ffi 



8^/6( g 2 +m2)( g 2 + m2) 



i w + (4m 2 + q\ + q\ + 2 qi • q 2 ) J 4(|q 1 + q 2 |) 



(B2) 



(B3) 



2. Two-pion one-pion-exchange 3N 



(Vi 



3N 



= 24 



' ,, s (kij -ki 3 ) 2 k i2 -ki 3 \ k i3,i?fu \ ( kl2 • kl3 ) 2 
-Fi(fci 2 ) F 2 (fci2j — — <r F 3 {k 12 ) — \-F4,{ki2) — 

-"13 ^13 -"13 -"13 



+ F 5 (k 12 )§^-F 6 (k 12 ) 



12)- 

ki2 • k 13 



K 13 u ^ "' K 13 
with structure functions F 1 (q) to F s (q) defined in Eqs. (2.17)-(2.20) of Ref. 



2F 7 (0)^+F 7 (k 12 )^L+AF s (k 12 )^-^ 

-"23 -"13 -"13 



(B4) 



3. Pion-ring 3N 



(V^ s ) = 8[9iJi(-ki2,ki 3 ) + 3fc 2 2 i? 2 (-k 12 , k 13 ) - 3k 12 • k u i? 3 (-k 12 , k 13 ) - 3k 12 • ki 3 i? 4 (-k 12 , k 13 ) 

+ 3fc 2 3 i? 5 (-ki 2 , k 13 ) - 6it6(-ki3, k 13 ) + 3i? 6 (-k 12 , k 13 ) - 2fc 2 3 i? 7 (-k 13 , k 13 ) + k\ 2 i? 7 (-k 12 , k 13 ) 
+ 2fc 2 3 i? 8 (-k 13 , k 13 ) - k 12 ■ k 13 i? 8 (-k 12 , k 13 ) + 2/c 2 3 i? 9 (-ki 3 , ki 3 ) - k 12 • k 13 # 9 (-k 12 , k 13 ) 
- 6i2io(-ki3, k 13 ) + 3i?io(-ki 2 , k 13 ) - 65i(-ki 2 , 0) + 3Sx(-k 12 ,k 13 ) + 3fc 2 2 S 2 (-k 12 , k 13 ) 
-3k i2 •ki 3 5 3 (-ki2,ki3) - 3k i2 • ki 3 S 4 (-ki 2 ,ki 3 ) + 3fc 2 3 5 , 5 (-ki 2 ,ki 3 ) + 95 6 (-ki 2 ,ki 3 ) 



with structure functions Ri and Si defined in Eqs. (A2) and (A7) of Ref. 



4. Two-pion-exchange contact 3N 
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Ct g A 
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5. Relativistic-corrections 3N 
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- (kia x k 13 ) ■ (ku x P 23 )F 1 3 /m (-k 1 2,k 13 ) - (ku x Pi 3 ) ■ (k 12 x k 13 )^ 1 4 /m (-k 12 ,k 13 ) 

- k 12 • ki 3 Ji B /m (-k 1 2,k 13 ,Pi2,P23,Pi3) + (kia x k 13 ) 2 J F 1 6 /m (-k 12 ,k 13 ) - k 12 • Pi 3 F 1 7 /m (-k 12 ,k 13 ) 



- ^2-Fl/m(-kl2,ki 3 ) + fcf 2 ^ 9 /m (-k 12 ,k 13 ) + /c 2 2 F 1 1 / ° m (-k 1 2,Pi2,P23) - k 12 • k 13 F^} m (k 12 ) 

-ki2 • k 13 F 1 1 / 2 m (fc 12 ) - k^ • m 23 F^f m (k 12 ) - k i2 • m l2 F^f m (k 12 ) 



(B7) 
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with 



Fu m (cii,<l2) 
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(qf + ml) 



(1 - 2&)(qi • q 2 ) 2 + (2~fo-l)qf 



^i 2 /m (qi>q2,q 3 ,q4) 



A 
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6. Three-pion-exchange and pion-interaction 4N 
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7. Two-pion-exchange contact 4N 
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